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INTRODUCTION 


This  report  concerns  the  analytical  studies  of  film  cooling 
in  gun  barrels  (Grant  No.  DAAG29“76*-G-0008,  Formerly  DAHC04--74-G--0191) 
during  the  period  October  1,  1974  to  December  31,  1976.  The 
studies  include  two  areas,  namely  cooling  with  wear  reducing  addi- 
tives  and  liquid  cooling.  The  dimensionless  wall  enthalpy'  dis¬ 
tribution  referred  to  as  film-cooling  effectiveness  has  been  obtained 
for  each  area.  In  this  section,  a  statement  of  the  problems  and 
a  list  of  publications  are  given.  In  the  ensuing  sections,  important 
results  from  the  analyses  will  be  summarized  followed  By  a  descrip-- 
tion  of  experimental  design.  Finally,  extension  of  the  present  work 
will  be  discussed. 

Two  areas  of  study  just  mentioned  are:  1)  Cooling  with  wear 
reducing  additives  (film  cooling  in  stationary  coordinates)  -  Consi¬ 
der  a  layer  of  the  ablative  polymeric  additive  wrapping  around  the 
propelling  charges.  When  the  high  temperature  gases  are  produced  during 
the  process  of  combustion,  the  additive  degrades  to  form  a  cool  gas 
film  that  flows  along  the  barrel  and  separates  the  high  temperature 
core  flow  from  the  barrel  surface.  The  cooling  effect  is  thus  ob¬ 
tained.  This  cooling  scheme  gives  the  high  film-cooling  effectiveness 
in  the  region  adjacent  to  the  position  of  shot  start  where  the  erosion 
is  most  series;  2)  Liquid  cooling  (film  cooling  in  a  moving  coordinate) 
In  this  cooling  scheme,  a  small  amount  of  liquid  is  confined  to, the 
inside  of  the  projectile.  Wlien  the  projectile  moves  under  high  pres¬ 
sure,  the  liquid  is  squeezed  out  to  cool  the  boundary  layer.  For  an 
observer  moving  with  the  projectile,  he  sees  that  the  coolant  flows 
out  from  a  fixed  location  into  the  unsteady  turbulent  boundary  layer. 

In' this  way,  the  problem  may  be  regarded  as  film  cooling  in  a  moving 


coordinate.  Liquid  cooling  technique  is  very  limited  in  the 
current  literature.  The  cooling  effect,  however,  appears  high. 

The  prediction  of  the  film-cooling  effectiveness  in  both,  two  areas 
is  our  essential  goal  for  this  project. 

There  are  three  publications  resulting  from  the  investi¬ 
gation  under  the  Grant.  They  are  listed  in  Table  I.  The  first 
two  and  third  in  the  table  treat  the  problems  in  the  first  area 
(cooling  with  wear  reducing  additives)  and  the  second  area  (liquid 
cooling)  respectively.  The  second  one  can  be  considered  an  improved 
version  of  the  first  one  since  the  former  gives  the  precise  treatment 


of  the  momentum  equation. 


I 


Table  I 

List  of  Publications 


1.  Tsou,  F.  K, ,  "Film  Cooling  in  Gun  Barrels",  presented  at  the  First 
International  Symposium  on  Ballistics,  Orlando,  Florida,  November,  1974 
and  published  in  the  proceedings, 

2.  Tsou,  F.  K. ,  "Prediction  of  Film-Cooling  Effectiveness  in  an  Unsteady 
Turbulent  Boundary  Layer,"  ASME  paper  No,  76-HT-28,  presented  at  the 
National  Heat  Transfer  Conference,  St.  Louis,  Mo.,  August,  1976. 

* 

3.  Tsou,  F.  K,  and  Varghese,  T. ,  "Film  Cooling  with  Injection  from  a  Moving 
Projectile,"  published  in  the  Proceedings  of  the  Second  International 
Symposium  on  Ballistics,  Daytona  Beach,  Florida,  March,  1976. 


Varghese  started  as  a  post-master  student  in  January,  1975  and  passed  his 
candidacy  examination  in  May,  1976.  He  is  currently  a  Ph.D,  candidate  in 
the  Department  of  Mechanical  Engineering  and  Mechanics,  Drexel  University. 
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Important  Results 


Prior  to  summarizing  the  results,  a  brief  description  of 

the  approach  used  in  the  investigation  is  given  here.  To  study 

the  problem  of  heat  transfer,  it  is  often  convenient  to  define 

a  heat  transfer  coefficient,  h  in  the  following  form; 

q  =  h  (T  -  T  .)  (1) 

To  utilize  this  definitive  expression,  one  needs  to  know  h  and 

the  reference  temperature,  If  the  adiabatic  wall  tempera^ 

ture,  T  is  used  as  T  h  for  the  cases  with  film  cooling  and 

without  film  cooling  is  approximately  the  same.  Since  h  for  the 

case  without  film  cooling  is  considered  available,  the  prediction 

of  T  thus  becomes  the  center  of  interest  in  a  film-cooling  pro- 
aw 

blem.  In  the  present  work,  the  dimensionless  enthalpy  known  as 
film-cooling  effectiveness  defined  below  is  used, 

::  ^  ~  (2) 

where  the  subscripts,  aw,",  i  are  referred  to  the  condition 
of  adiabatic  wall,  free  stream,  and  injection  respectively.  The 
determination  of  n  is  therefore  the  essential  undertaking  of  this 
project.  Employing  typical  assumptions,  the  governing  integral 
equation  can  be  recast  into  three  first  order  equations  that  are 
hyperbolic.  The  method  of  characteristic  is  then  applied  to  ob¬ 
tain  the  solution  of  p.  The  results  of  calculation  from  the  publi¬ 
cations  mentioned  in  the  foregoing  are  as  follows : 

1)  Cooling  with  wear  reducing  additives 

Fig.  1  shows  the  effectiveness  p  versus  the  dimensionless 
distance  with  time  in  milli-second  as  a  parameter.  A3.1  curves 

start  from  the  position  of  shot  start,  At  early  times  (t  =0.2 

0.3  msec),  n  decreases  monotonically .  After  0.3  msec,  the  rate"  of 
ablation  decreases  and  reduces  to  zero  at  t  =  0.5  msec.  Thus  for 

t  =  0.4  msec,  increases  p  first  and  then  decreases.  The  curves  for 
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t  =0.5  and  0.7  msec  have  the  similar  nature.  At  t  =  0,5  msec, 

Ti~  0  at  =  1*  The  effectiveness  n  at  t  =  0.7  msec  remains 

zero  until  x/x^  =  1.4.  These  results  indicate  that,  at  earlier 
times,  the  cooling  effect  is  very  significant  in  the  region 
where  the  erosion  is  most  series. 

2)  Liquid  Cooling 

A  similar  plot  in  a  stationary  coordinate  is  shown  in  Fig.  2. 

The  curves  again  were  plotted  from  the  position  of  shot  start.  At 
t  =  0,2  and  0,5  msec,n  increases  and  reaches  to  unity  at  the 
base  of  the  projectile.  The  third  curve  at  t  =1  -msec  shows  a  gradual 
decrease  first  and  then  an  increase  to  -11=1  at  the  base.  Like 
curves  in  Fig.  1,  there  is  a  point  of  discontinuity  on  each  curve. 

The  gradual  decrease  in  the  region  of  small  x  is  probably  caused  by 
the  assumptions  used  in  the  analysis.  Generally  speaking,  the  effec¬ 
tiveness  of  liquid  cooling  is  fairly  high.  However,  in  the  problems 
of  gun  barrel  cooling,  one  requires  better  cooling  effect  near  the 
position  of  the  shot  start.  The  trend  shown  in  Fig.  2  is  just  op¬ 
posite. 

With  the  results  of  the  effectiveness  just  mentioned,  one  is 
able  to  calculate  the  erosion  of  a  surface.  For  example,  in  the  case 
of  one-dimensional  melting  of  a  solid  with  complete  removal  of  meet, 
Goodman’s  integral  method  can  be  conveniently  applied.  The  required 
surface  heat  flux  in  this  method  is  obtained  from  Eq.  (1),  in  which, 
the  heat  transfer  coefficient  h  is  given  and  the  adiabatic  wall  tempera¬ 
ture,  T  is  computed  from  the  present  solution. 
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Experimental  Design 


Experimental  work  using  gun  barrels  are  available  in  the 
literature.  The  fact  that,  in  the  gun  barrel  measurements,  the 
ablative  rate  of  the  additives  is  difficult  to  control  preclude 
comparison  v/ith  the  present  analytical  results.  To  substantiate 
these  results,  it  is  felt  that  the  use  of  an  intermittent  subsonic  wind 
tunnel,  known  as  Ludweig  Tube,  would  be  appropriate.  The  tube  con¬ 
sists  of  a  conventional  shock  tube  with  the  insertion  of  a  converging- 
diverging  nozzle  into  the  section  upstream  of  its  diaphragm.  A 
transient  subsonic  flow  in  a  period  of  several  milli-seconds  can  be 
obtained  in  the  test  section.  The  injection  system  including  a  quick 
closing  valve  mechanism  will  give  desired  injection  mass  flow.  In 
this  way,  the  transient  film-cooling  effectiveness  can  be  measured. 

Some  details  of  the  experimental  method  are  given  in  Appendix.  The 
method  applied  to  the  condition  that  is  relevant  to  cooling  with  x^^ear 
reducing  additives . 

Discussion  and  Conclusion 

Extension  of  the  present  work  other  than  the  experimental  measure¬ 
ments  just-described  is  given  here.  The  discussion  is  again  limited 
to  the  case  of  cooling  with  wear  reducing  additives. 

The  problems  that  remain  to  be  studied  are:  1)  The  boundary 
la^^er  may  be  destroyed  initiall^^  so  as  to  give  no  cooling  effect. 

2)  To  study  the  ablative  rate,  one  has  to  consider  the  parameters: 
the  location,  length,  and  thickness  of  the  linear.  Some  nev7  models 
or  modification  of  the  present  models  need  to  be  developed.  3)  The 
density  of  the  ablative  materials  may  be  another  parameter.  .A  study 
based  on  the  diffusion  thermal  effect  would  give  an  ansx<rer.  4)  The 
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use  of  flaps  has  shown  in  improvement  in  cooling  effect. 

Analytical  study  of  transient  flow  and  heat  transfer  over  such 
flaps  can  be  made. 

The  ultimate  goal  for  the  present  studies  is  to  obtain 
information  on  the  inert  erosion  of  gun  barrel  surfaces.  As  men¬ 
tioned  earlier,  it  is  convenient  to  solve  one  dimensional  phase 
change  problems  with  a  proper  numerical  scheme.  Since  the  axial 
heat  conduction  is  significant,  it  would  be  more  important  to  solve 
two-dimensional  problems. 

Tlie  approach  used  in  the  present  study  may  have  other  appli¬ 
cations,  e.g.  film  cooling  along  gas  turbine  blades.  Here,  the 
problems  of  interest  include  the  process  of  start-^up,  shut-down 
and  the  period  of  load  change  as  well  as  the  wake  effect.  The 
transient  response  of  blade  surface  temperature  will  contribute 
to  the  theriii.al  fatigue  and  creep  of  the  blade  material  significantly. 


Appendix  f 


This  appendix:  describes  hew  expeiriinental  de  terrain  at  ion  of 
film-” cooling  effectiveness  can  be  perfomiecL  Consider  gas  injec¬ 
tion  into  the  subsonic  turbulent  boundary  layer  on  a  flat,  adia¬ 
batic  vjall.  The  following  six  combinations  of  injection  gas  mass 
flow  rate,  and  free  stream  gas  velocity,  can  be  studied. 


(1) 

0 

m .  = 

const . 

11 

const . 

(2) 

0 

m .  = 

1 

m^(t) 

U  = 

oo 

const . 

(3) 

0 

m .  = 

1 

const . 

u  - 

CO 

(4) 

0 

m .  = 

1 

Sgt) 

u  = 

CO 

(5) 

0 

m .  = 

1 

const . 

u  = 

oo 

(6) 

0 

m. 

X 

\(t) 

u  = 

CO 

U^(x,t) 

,  Each  of  these  cases  is  qualitatively  shown  in  Fig.  Al.  The  tech¬ 
niques  used  to  obtain  free  stream  flow  and  injection  flow  as  well 
as  film-cooling  effectiveness  measurements  are  described  in  the 
following : 

Subsonic  flow  production  technique.  Use  is  made  of  an 
intermittent  subsonic  vrind  tunnel  that  consists  of  a  conventional 
shocl^  tube  with  the  insertion  of  a  converging-diverging  nozzle  into 
the  section  upstream  of  the  diaphragm  (high  pressure  section) .  Such 
.  a  device  is  knov-n  as  Ludweig  Tube.  The  operating  stages  of  this  tube 
are  shov/n  in  Fig.  A2 ,  an  x-t  diagram  and  the  sketch  of  the  tube. 

Referring  to  Fig.  A2 ,  vfnen  the  diaphragm  is  broken  at  time 
zero,  a  shock  wave  follovjed  by  a  contact  surface’  travels  dov.mstreaiu 
into  the  low  pressure  gas.  Upstream  of  the  diaphragm,  an  expansion 

t  The  author  is  indebted  to  Professor  L.  T.  Smith,  Department  of 
Mechanical  Engineering  and  Mechanics,  Drexel  University,  for  his^ 
important  contribution  to  this  appendix. 
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wave  propagate  into  a  high  pressure  gas,  accelerating  it  toward 
the  nozzle  exit.  During  this  initial  starting  period,  we  have 
Uoo  *=  Uoo(x,t)  in  the  subsonic  test  section  of  the  nozzle,  and  at 
a  given  location,  the  variation  of  velocity  with  time  is  as  shown 
in  Fig.  A1  (a  and  b) .  As  soon  as  the  nozzle  throat  becomes  choked, 
as  determined  by  the  ratio  of  nozzle  throat  area  to  '^supply 
tube  area,  the  remaining  part  of  the  expansion  wave 
is  swept  downstream.  In  the  supply  tube,  the  strength,  or  pressure 
ratio,  of  the  expansion  wave  becomes  fixed.  Thus,  a  region  of  con¬ 
stant  gas  properties  and  velocity  is  created  between  the  expansion 

^  I 

wave  and  the  nozzle  inlet.  This  region  serves  as  the  supply  con¬ 
dition  for  a  steady  expansion  flow  through  the  nozzle,  with  Ho  = 
const,  in  the  subsonic  test  section  as  shown  in  Fig.  A1  Cc  and  d)  . 
Steady  flow  lasts  until  the  expansion  wave  reflectes  from  the  end 
wall  of  the  supply  tube  and  returns  to  the  nozzle  inlet.  At  this 
stage  in  the  operation,  the  expansion  wave  reflects  from  the  nozzle 
and  propagates  back  upstream  into  the  supply  tube  to  produce  a  sec¬ 
ond  period  of  steady  nozzle  flow  with  altered  free  stream  conditions 
Subsequent  expansion  wave  reflections  between  tthe  nozzle  and  supply 
tube  end  wall  produce  additional  periods  of  steady  nozzle  flow  as  11 
lustrated  in  Fig.  A1  (e  and  f )  .  We  thus  have  the  case  (t)  . 

Finally,  a  reflected  shock  system  returns  to  the  nozzle  from  the  low 
pressure  tube  to  destroy  the  periodic  nozzle  flow  pattern. 

For  each  of  these  experiments,  nozzle  freie  stream  velocity 
and  gas  properties  will  be  computed  using  existing  theory.  These 
calculations  will  be  verified  by  comparing  measured  static  pressure 
in  the  test  section  to  predicted  values.  Mori&over,  measurements  of 
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the  static  pressure  will  be  used  to  adjust  the  theoretical 
calculations  to  account  for  the  effects  of  the  viscous  boundary 
layer  and  thus  arrive  at  an  accurate  determination  of  the  free 
stream  conditions  in  the  test  section. 

For  the  initial  transient  flow  period,  with  t), 

the  governing  equations  of  motion  will  be  solved  numerically, 
using  the  method  of  characteristics,  to  determine  free  stream 
conditions  as  a  function  of  time  at  a  given  location  in  the 
nozzle  test  section  11,2].  Existing  analytic  solutions  13]  will 
be  used  to  compute  test  times  and  free  stream  conditions  during 
the  initial  steady  flow  period,  =  const.,  and  during  the  peri¬ 
odic  flow,  U„  =  (t)  . 

C2)  Gas  Inlection  Technique.  The  gas  injection  system  can 
be  modeled  after  the  Oxford  system  14,  5],  however,  certain  modi¬ 
fications  have  to  be  made  in  order  to  achieve  transient  injection 
with  m°  =  m^(t).  Referring  to  Rig.  A3,  injection  gas  will  be  held 
in  a  chamber  of  volume  V  at  pressure  p^.  Before  the  start  of  an 
experiment,  the  chamber  will  be  sealed  off  from  the  channel  leading 
to  the  injection  slot  in  the  nozzle  test  section  by  means  of  a  thin 
Mylar  diaphragm.  Injection  will  be  started  by  bursting  the  diaphragm  with  a 
needle  activated  by  a  solenoid.  A  current  pulse  of  3000  amps .with  a  rise  time 
200  ps  [4,5]  will  be  used  to  activate  the  solenoid,  and  will  be  trig¬ 
gered  by  a  pressure  pulse  generated  by  the  shock  wave  which  travels 
into  the  low  pressure  dump  tube  section  of  the  tunnel.  A  time  delay 
system  will  be  provided  to  control  the  time  of  the  start  of  gas  in¬ 
jection,  and  the  design  will  allow  injection  to  be  started  prior  to 
the  free  stream  flow  in  the  test  section. 
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For  steady  gas  injection,  =  const.,  the  injection  system 
will  be  operated  without  the  quick  shut-off  valve  shown  in  Fig.  A3. 

The  injection  mass  flow  rate  will  be  determined  by  measuring 
the  time  rate  of  change  of  the  injection  pressure  dp^/dt,  during 
the  run  [4].  For  the  short  injection  times  involved  (10  to  15  milli¬ 
seconds)  ,  the  process  within  the  injection  gas  reservoir  can  be  as¬ 
sumed  adiabatic.  The  injection  mass  flow  rate  is  then  given  by 

dP. 

0  _  y _  _ X 

“i  yRT^  dt 

where  y  is  the  ratio  of  specific  heats  of  the  injection  gas,  R 
Is  the  gas  constant,  and  is  the  injected  gas  temperature.  For 
small  reservoir  pressure  drops,  the  injection  gas  temperature,  T^, 
may  be  considered  constant  14].  The  gas  temperature  at  the  point 
of  injection  into  the  turbulent  boundary  layer  is  taken  to  be  equal 
to  T^,  and  this  assumption  can  be  checked  using  a  thin  film  tempera¬ 
ture  gage  mounted  flush  to  the  channel  conducting  gas  from  the  reser¬ 
voir  to  the  injection  slot,  as  shown  in  Fig.  A3. 

For  transient  gas  injection,  m^  =  m^(t) ,  a  quick  shut-off 
valve  will  be  mounted  on  the  diaphragm  bursting  needle  as  shown  in 
Fig.  A3.  With  this  configuration,  one  can  first  have  several  milli¬ 
seconds  of  steady  injection  followed  by  a  rapid  decrease  in  injection 
mass  flow  rate  with  time  as  shown  in  Fig.  A1  (b,  d  and  f ) . 

(3)  Film  cooling  effectiveness  measurements.  During  the  short 
testing  times  involved  in  the  operation  of  the  Ludwieg  tube,  there 
is  not  enough  time  for  appreciable  heat  transfer  into  the  test  section 
nozzle  blocks.  We  thus  deal  x^ith  the  case  of  injection  into  the  tur- 
bulend  boundary  layer  over  an  adiabatic  wall.  The  effectiveness  has 
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been  defined  as  the  dimensionless  adiabatic  wall  enthalpy  (see  Eq.  2 
in  the  text).  For  the  case  of  injection  gas  and  free  stream  gas  of 
the  same  chemical  species,  and  assuming  both  to  be  an  ideal  gas  with 
constant  specific  heats,  the  effectiveness  becomes 

T  -  T 

,  s  aw  “ 
ri(x,t)  =  — 17^ — 

i  “ 

where  T  is  the  adiabatic  wall  temperature,  and  and  T  are  the 
aw 

temperatures  of  the  free  stream  gas  flow  and  injection  gas  respectively. 

With  both  T  and  T  determined  by  the  methods  discussed  previously,  it 
i  “> 

remains  to  experimentally  determine  adiabatic  wall  temperature  in. order 
to  determine  n. 

Platinum  thin  resistance  gages  will  be  used  to  measure  adiabatic 
wall  temperature  as  a  function  of  time  at  several  locations  upstream  and 
downstream  of  the  gas  injection  port.  The  gages  16]  consist  of  a  pyrex 
or  quartz  insulating  material  mounted  flush  with  the  wall  surface.  The 
surface  temperature  of  the  Insulating  material  is  measured  by  a  thin 
platinum  film  in  good. thermal  contact  with  the  surface.  The  thickness 
of  the  film  is  of  the  order  of  one  micron.  A  constant  current  is  passed 
through  the  film  such  that  measurement  of  the  voltage  across  the  film 
enables  the  resistance  and  hence  the  surface  temperature  to  be  deter¬ 
mined.  The  rise  time  of  such  gages  are  typically  of  the  order  of  10 
microseconds  [4] ,  two  orders  of  magnitude  less  than  the  expected  testing 
time  of  the  experiments . 
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Figure  Al.  Combinations  of  free  stream  velocity,  U^,  and  injection  gas  mass 


flow  rate,  m. ,  versus  tine.  Typical  time  scales  shown  are: 
=  4  ms,  t2  =  16  ms,  t^  =  70  ms. 
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Figure  A2.  The  Ludwieg  Tube  intermittent  wind  tunnel:  sketch  and  simplified  x-t  diagram 
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